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Abstract 

We calculate the masses of the ++ , 0~~ and 1 H glueballs from 
QCD in 3+1 dimensions using an eigenvalue equation method for 
Hamiltonian lattice QCD developed and described in j|, ||, |(| by the au- 
thors. The mass ratios become approximately constants in the coupling 
region 6/g 2 € [6.0,6.4], from which we estimate M(0 — )/M(0 ++ ) = 
2.44 ±0.05 ±0.20 and M (l+~)/M (0++) = 1.91 ± 0.05 ± 0.12. 

PACS number(s): ll.15.Ha, 12.38.Gc, 14.20.-c,14.40.-n 
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1 INTRODUCTION 



If QCD is the correct fundamental theory of strong interactions, a 
lot of glueball states should exist. These glueballs are gluonic bound 
states formed by strong self-interactions of gluons, and characterized 
by various quantum numbers J PC . Precise calculation of the glueball 
masses will distinguish QCD from other models and give important 
guide for experimentalists. Among various nonperturbative techniques, 
lattice QCD seems to give more reliable estimates since it is model- 
independent. 

Within lattice QCD, Monte Carlo simulation on an Euclidean lattice 
is still the most popular method, since in the Euclidean formulation, 
space-time symmetry is apparent, the particle spectrum is given by the 
correlation length of Green's function, and one can use the simulation 
methods in statistical mechanics. In fact, it can easily bridge the strong 
coupling region where confinement is obvious, and weak coupling region 
which is more physically relevant. With increasing lattice sizes, higher 
statistics and more suitable techniques, it might evidently give reliable 
predication for the glueball masses. For a review of the current status, 
see ref. |TJ. 

The Hamiltonian formulation is an alternative to the standard Eu- 
clidean approach. In such a formulation, one works with a conventional 
quantum mechanical system, and the spectrum and wavefunction cor- 
respond to an eigenvalue and eigenstate of the Hamiltonian. Further- 
more, it might possibly give new insight into the structure of the the- 
ory. Over the last twenty years, many attempts have been made. The 
main difficulty of some conventional methods (e.g. strong coupling 
expansion) is that they converge very slowly and very higher order 
1/g 2 calculations are required to extend the results to the intermediate 
coupling region. Unfortunately, high order calculations are difficult in 
practice. For some other recent and interesting developments, see 0, |3[] 
for details. 

In refs. [|, [5], |6|], we proposed an efficient eigenvalue equation 
method for Hamiltonian lattice gauge theory. This method consists of 
solving order by order the eigenvalue equation with a suitable scheme 
preserving the correct continuum behavior at any order of approxi- 
mation. So far, we have applied this method to U(l), SU(2) and 
SU(3) models in 2+1 dimensions and nonlinear O(N) o model in 2 
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dimensions. The reliability of our method is supported by increasing 
evidences @, § §, 0, g, |, 0, 



Most of these results have been 



summarized in [|T2], pf ]. 

In a recent paper 
QCD 3 , more accurate than the previous papers 



we obtained the lightest glueball mass in 

with finite order 



truncation error under well control, 
sional reduction 



m 16 



We then used the idea of dimen- 
to estimate the ++ glueball mass in four 
dimensional QCD. Our results suggested M(0 ++ ) ps 1.71 ± 0.05 Gev, 
in nice agreement with the most recent Monte Carlo data M(0 ++ ) ~ 
1.740 ± 0.071 Gev by the IBM group |7|. This favors 0//j(171O) as a 
candidate of the ++ glueball. 

It is very desirable to do concrete computations in the realistic 
theory: QCD in 3+1 dimensions. In this paper, we try to make a 
first step towards this direction. The rest of the paper is organized as 
follows. In Sec. [|, we recapitulate the basic elements of our method. 
The results for the glueball masses are given in Sec. |3[ In the Appendix. 
0, we discuss the relation between the Hamiltonian and Euclidean 
formulations for the purpose of comparison. 



2 OUR METHOD 



Our method has been described in detail in refs. ||, |5|, ^ |7|, [^, [| [II], O, 
12], [13]. For the reader's convenience, we will repeat the most important 



points. The Schrodinger equation H\Q) = 6q\Q) on the Hamiltonian 
lattice for the ground state 

\n)=exp[R(U)]\0) (1) 

and vacuum energy can be reformulated as 

[E h R{U)}} + [E h R{U)][E h R{U)}} 



-^Y,Tr(U p + Ui) = 2 ^en. (2) 

y p y 

To solve this equation, let us express R{U) in order of independent 
set of gauge invariant operators G n ^, i.e., R{U) = J2 n Rn(U) = 
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J2 n ,i C n ,iGn t i(U). Substituting it to (0), we have the iVth order trun- 
cated eigenvalue equation ||. |5[ 

N 

J2{[Ei,[Ei,J2MU)}}+ E [E h R ni (U)][E h R n2 (U)]} 

l n m+n 2 <N 

-±J2ReTr(U p ) = ^e u . (3) 
9 p 9 

Since the higher order graphs are created by the second term in ([|), 
at Nth order approximation all contributions from this term with n\ + 
112 > N should be discarded. The reason for such a truncation method 
is to maintain the correct continuum behavior during the truncation, as 
explained in detail in ref. ||]. By taking the coefficients of the graphs 
G n> i in this equation to zero, we obtain a set of non-linear algebraic 
equations, from which C n ^ are determined. Therefore, solving lattice 
field theory is reduced to solving the algebraic equations. Similarly, the 
eigenvalue equation H\F) = ep\F) for the glueball mass M(J PC ) = 
Ae = ep — £q and its wave function 

,^)-<ra],,> 

at order iV is 

N 

J2{[Ei,[E l ,Y,F n (U)}}+2 £ [E h F ni (U)}[E h R n2 (U)]} 

l n m+n 2 <N 

= ^F.(V), (5) 

9 n 

from which the coefficients f n ^ in the glueball operator F(U) = 
J2 n E n (U) = J2 n ,ifn,iGni are determined. According to our experi- 
ences, it is hopeful that the results converge rapidly to some stable 
value. The difference between the results from different truncation 
orders is our estimate for the systematic error. 

Other detailed techniques, like the choice of independent set of op- 
erators in SU(3) and efficient ways to improve the continuum behavior 
of the physical quantities, can be found in refs. p|, p|, [TT . 
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3 RESULTS AND DISCUSSION 



Since this work is our first attempt to apply the method above to 
3+1 dimensional pure QCD, we would like to compute the masses of 
glueballs with gluonic operators easily constructed. Here we choose 
three glueballs ++ , , and Some lowest order gluonic operators 
are plotted in Fig. 1. The operators R n and F n are linear combinations 
of these graphs, and the first order ones are 

R^diG^ + G^ + G^ + h.c, 

Ff- =fl + -(Gf\-h.c). (6) 

Since the ++ glueball has symmetry the same as the vacuum, we take 
F® ++ = R\. Once the lowest order graphs are determined, the higher 
order graphs are generated by the second term in (^) for R n or that 
in (|^) for F n , from which an independent set is chosen, taking into 
account the uni-modular conditions. The operators so generated are 
certainly of the symmetry of the corresponding glueballs. For example, 
at the 2nd order 

P 2 — J2,l 1^2,1 + ^2,1 + ^2,1) + J2,2 1^2,2 + ^2,2 + ^2,2^ 
+/2,3 1^2,3 + ^2,3 + ^2,3^ 
+ i2,4 1^2,4 + ^2,4 + ^2,4^ + J2,5 1^2,5 + ^2,5 + U 2,5) 

+J2.6 l (jr 2,6 + ^2,6 + ^2,6) 

, fO++ (r (l) , r (2) , ^(3) ~(4) ^(5) ^(6) , ^(7) ^(8) ^(9) 
'J2,7 \ U 2J i ( - T 2,7 + ^2,7 + ^2,7 + < - jr 2,7 + ( - T 2,7 + ( - T 2,7 + ^2,7 + < ~ T 2,7 

,^(10) , ^(H) , r (12), 
+ Lr 2)7 + (_r 2 ,7 + Cr 2 ,7 J 

, f 0++ (r (l) , ^(2) ^(3) ^(4) ~(5) ~(6) ^(7) ^(8) ^(9) 
+ i2,8 1^2,8 + ^2,8 + ^2,8 + ^2,8 + ^2,8 + ^2,8 + ^2,8 + ^2,8 + ^2,8 
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+ Lr 2 ,8 + ^2,8 + ^2,8 ) 

, f 0++/^(l) , r (2) , r (3) , r (4) , r (5) , r (6) , r (7) , r (8) , r (9) , ,,(10) 
+ i2,9 l"2,9 + ^2,9 + ^2,9 + ^2,9 + ^2,9 + ^2,9 + ^2,9 + ^2,9 + ^2,9 + ^2,9 

+Lr 2 ,9 + ^2,9 ) 

, f 0++/Wl) ,/^(2) 1/^(3) , r (4) , r (5) r (6) 
"rJ2,10 l"2,10 ' "2,10 ~t~ "2,10 t "2,10 ~r "2,10 ~t~ "2,10 

,^(7) , r (8) , r (9) ^(10) ^(11) r (12h , 
+ ^2,10 + "2,10 + "2,10 + "2,10 + "2,10 + "2,10j + n - C - 

p0~ _ f 0-/^(1) , ^(2) , ^(3) ^(4) ^(5) ^(6) 
-^2 — i2,7 l>"2,7 + ^2,7 + ^2,7 + ^2,7 + ^2,7 + ^2,7 

+ Lr 2 ,7 + ^2,7 + ^2,7 + ^2,7 + ^2,7 + ^2,7 ) 

, fO— fr (l) r (2) r (3) r (4) r (5) r (6) r (7) r (8) r (9) 
+J2,8 \ U 2,8 + ^2,8 + ^2,8 + ^2,8 + ^2,8 + ^2,8 + ^2,8 + ^2,8 + ^2,8 

, ^(10) , Mil) M 12 h 

+J2,9 l"2,9 + ^2,9 + ^2,9 + ^2,9 + ^2,9 + ^2,9 + ^2,9 + ^2,9 + ^2,9 

+ Lr 2 9 + Lt 2 9 + Lr 2 ,9 ) 

, fO—fM 1 ) ±r( 4 ) _i_r'( 6 ) 

"TJ2,10 l"2,10 "r "2,10 "r "2,10 "r "2,10 "2,10 "r "2,10 

, r (7) , r (8) ^(9) ^(10) ^(11) r (12)x , 
+ ^2,10 + "2,10 + "2,10 + "2,10 + "2,10 + Lt 2,1oJ — n - C - 

_ fl+-/-<(3) , f l+"^(3) , f l+- r (3) 
^2 — /2,1 tjr 2,l + J2,3 U 2,3 + /2,4 U 2,4 



+■72,7 l tjr 2, 7 + < - r 2,7 + ( - J 2,7 + ( - T 2,7 + ^2,7 + ^2,7 ' ( - T 2,7 + ( - T 2,7 ) 

, f 1+- ^(5) , r (6) ^(7) ^(8) ^(9) ^(10) r (ll) r (12)x 
+J2,8 1^2,8 + Lr 2,8 + Lr 2,8 + Lr 2,8 + ^2,8 + ^2,8 + Lr 2,8 + ^2,8 J 

, fl+~ ( r (JS) , ^(6) , ^(7) ~(B) , ^(9) ^(10) r (ll) r (12)x 
+i2,9 l (jr 2,9 + ^2,9 + Lr 2,9 + ^2,9 + ^2,9 + ^2,9 + ^2,9 + ^2,9 ) 

>fi+-(r® 4_r( 7 ) -u^ 8 ) , r( 10 ) , r^-L^h 

"TJ2,10 1 Lt 2,10 i (jr 2, 10 i* ^2, 10 "T (jr 2, 10 + ^2,10 + ^2,10 ' ^2,10 + Lt 2,10J 

- /i.e. (7) 

Higher order graphs can be calculated in a similar way. 

At relatively strong coupling, the absolute value of oMg, calculated 
in the Hamiltonian formulation and converted to the Euclidean one 
using (|IBD, differs from the result from strong coupling expansion in 
the Euclidean formulation. This is not surprising because the Hamilto- 
nian and Euclidean formulations are different schemes at finite lattice 
spacing, and the weak coupling relation (|l^) doesn't not hold in the 
strong coupling region. For (3 > 6.0, we do observe that the differ- 
ence in aM.E becomes much smaller at weaker coupling and there is a 
tendency approaching the Monte Carlo data 0. 

Similar to the most recent Monte Carlo data [|l[] in the available 
coupling region, clear asymptotic scaling window for the individual 
mass aM(0 ++ ), aM(0 ), or aM(l + ~) could not be found. There 
might be two possible reasons for this scaling violation: the available 
coupling region being not weak enough for the asymptotic scaling re- 
lation between a and (3e to be valid or the results being not accurate 
enough. The first one may be reduced by the Symazik's improvement 
or Lepage-Mackenzie scheme. The second one may also be improved 
by higher order calculations. For the mass ratio Mjp/Mjp, the large 
part of errors in M E 1] and il4 2) are canceled. Indeed, we observe ap- 
proximate plateau for the mass ratio both for M(0 )/M(0 ++ ) and 
M(l +_ )/M(0 ++ ) in the coupling region (3 G [6.0, 6.4], which are shown 
in Fig. 2. 

From the plateaus in (3 G [6.0, 6.4], we estimate 
Af(0~) 

— V^tt = 2.44 ± 0.05 ± 0.20 
M 0++ 
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(8) 



where the mean value is the averaged one over the data in this region, 
the first error is the error of the data in the plateau, and the second 
error is the upper limit of the systematic errors due to the finite order 
truncation. 

For comparison, we list corresponding results from other lattice cal- 
culations: For M(0 )/M(0 ++ ), it is approximately 3.43 ± 1.50 from 
Monte Carlo simulations (see the references in [l|]), 2.2 ± 0.2 from t- 
expansion plus Pade approximants [19]. For M(l + ~)/M(0 ++ ), it is 
1.88 ± 0.02 from Monte Carlo, and 1.60 ± 0.60 from t-expansion plus 
Pade. For M(0 )/M(0 ++ ), we obtain a value between the Monte 
Carlo and t-expansion data, with error under better control than the 
former one. For the M(1 H )/Af(0 ++ ), where the systematic error in 
the Monte Carlo data is very small, we get a value consistent well with 
the Monte Carlo data. 

In conclusion, we have done spectrum computations of the ++ , 
and 1 H glueballs, by means of the eigenvalue equation method in 
& |5|, ||. Our first results for 3+1 dimensional QCD from this method 
are encouraging. The spectrum calculations of other glueballs, such as 
2 ++ , 0~ + and + ~, are in progress. 
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A EUCLIDEAN AND HAMILTONIAN FOR- 
MULATIONS 



The standard Yang-Mills action is 

S = ~\j d A xF« u F« u . (9) 
On the Euclidean lattice, the standard Wilson action is 

S = -^Y,Tr(U p + Vl-2), (10) 

where (3% is related to the number of color N c and the bare coupling 
9e by Pe = 2N C / g%, and U p = is the product of four link variables 
around an elementary plaquette in the /i and v directions. In the 
continuum limit (3e — » oo, it reduces to the Yang-Mills action (§). 

To derive the lattice Hamiltonian from (fL0|) , let us fix the lattice 
spacing a s = a in spatial directions and decrease the lattice spacing a t 
in the temporal direction so that £ = a/a t — > oo. In this case, (|T0|) 
becomes 

S = - A £ Tr(U u + h.c. - 2) - A £ Tr(C/ . . + fc. c . _ 2 ). (11) 

c i c i>j 

Quantum- mechanically, 

2iV c 2iV c 
St 9 2 s 

where g t ^ g s due to the £ factor and quantum effects. They are related 
to the Euclidean coupling by 

\ = ^ + c t + 0(g 2 E ), \ = \ + c s + 0{g E ) (13) 
91 9e 9 s 9e 



if gE is small enough. The coefficients are calculated in fig] : 

1 W173 

ct = 4iV c [-0.01631 + — ], c s = 4iV c [0.01707 - -^p]- (14) 
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Using the transfer matrix or Legendre transformation methods, one 
can show that the lattice Hamiltonian in the temporal gauge is 

H = ^H KS , (15) 
9s 

where 

Hks 4E«) - -\j2 Tr (U P + Ul -2). (16) 
za x ,j a 9u v 

is the Kogut-Susskind Hamiltonian and 

9h = V9t9~ s (17) 

is the gauge coupling of the Kogut-Susskind Hamiltonian gauge theory. 
In (|16D, the summation J2x,j is over the link / in the positive spatial 
directions, and the summation J2 P is over the spatial elementary pla- 
quette. For weak enough coupling qe — > 0, the mass can be compared 
with the Euclidean one 

M E = Mh— (18) 
9s 

if gu is converted to g& by (|17D . The gt/g s factor doesn't matter if one 
calculates the mass ratio. 
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Figure Captions 

Fig.l Some relevant gluonic operators at the lowest orders. 

Fig. 2 Mass ratios of the glueballs. 

12 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9609435vl 



This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9609435vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9609435vl 



This figure "fig2-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9609435vl 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9609435vl 



